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Performance improvement was attained in data reconstructions of 2-dimensional tunable diode laser absorption 


spectroscopy (TDLAS). Multiplicative Algebraic Reconstruction Technique (MART) algorithm was adopted for 


data reconstruction. The data obtained in an experiment for the measurement of temperature and concentration 


fields of gas flows were used. The measurement theory is based upon the Beer-Lambert law, and the measurement 


system consists of a tunable laser, collimators, detectors, and an analyzer. Methane was used as a fuel for com- 


bustion with air in the Bunsen-type burner. The data used for the reconstruction are from the optical signals of 


8-laser beams passed on a cross-section of the methane flame. The performances of MART algorithm in data re- 


construction were validated and compared with those obtained by Algebraic Reconstruction Technique (ART) al- 


gorithm. 
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Introduction 


Temperature distribution on the cross-section of a 
combustion flame allows for exact analyses on the com- 
bustion phenomena. Measurements on the exhaust gases 
such as NOx, SOx, COx and CHx from plants or facto- 
ries allows for preventing air pollutions and global 
warming. TDLAS (tunable diode laser absorption spec- 
troscopy) has been used for these purposes. 

For the measurement on the concentrations of the gas- 
es, there have been several representative pointwise tech- 
niques. Wang et al. [1] used a WAV (surface acoustic 
wave) type sensor to detect CO, and humidity. Stemha- 
gen et al.[2] used a WAV sensor to detect gas and tem- 
perature. Fine et al. [3] constructed a semiconductor-type 
sensor to detect CO, NOx, NH3 and CO, gases. Chen [4] 
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proposed a thermal conductivity-type sensor for detection 
of H; gas. Dossi et al.[5] adopted an electrochemical-type 
sensor for NH; gas. Since all of these conventional mea- 
surement techniques were based on pointwise measure- 
ments, the distribution of temperature or concentration 
can’t be measured. As field-wise measurement tech- 
niques, Ko et al. [6] proposed a digital speckle method 
which allows for the measurements of spatial distribution 
of CO, concentration. Since gas concentration is affected 
by the gas temperature, temperature distribution should 
be measured to predict the gas concentrations. Deguchi et 
al. [7] proposed a technique which can quantify the con- 
centration distributions of HO vapor by measuring the 
temperature distribution of a gas flame. They used a tun- 
able laser and introduced 2-dimensional array of the laser 
beam, in which 8 laser beam lines were installed and the 
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intensity of the laser lines were sensed by 8 detectors 
(photodiode). The intensity signals from 8 detectors were 
used to reconstruct the spectrum patterns of a gas. They 
used the ART (algebraic reconstruction technique) [8] 
method for reconstructing the concentration and temper- 
ature fields. In the meantime, Doh et al. [9] developed a 
tomographic PTV (particle tracking velocimetry) in 
which spatial distribution of velocity vectors are obtained 
from a set of tomography images that were reconstructed 
by the MART (multiplicative algebraic reconstruction 
technique) method. In this study, it was known that the 
MART algorithm allows less calculation load and less 
calculation errors than the ART algorithm. In this study, 
performance of the MART algorithm is shown by using 
the experimental data obtained by the Deguchi group [7]. 


Temperature vs Relative Intensity 


Fig. 1 shows major spectroscopic steps, such as ab- 
sorption, fluorescence and photoionization, from which 
spectra signals can be measured for the investigations of 
gas properties. The absorption spectroscopy is funda- 
mentally based upon Beer-Lambert Law as shown in Fig. 
2. When the incident laser beam passes through the uni- 
form medium, the ratio of the transmitted light intensity 
(I) against the incident light intensity (Zo) is 

I,(A) = [p(A)-exp{-@} (1) 
inversely proportional to the power of the absorbance 
coefficient (æ, ) Eq. (1). The intensity ratio of incident 


light and transmitted light depends on the mole fraction. 
The number density of the measured species n is related 
to the amount of light absorbed as in the following for- 
mula. 

L(A)/[p(A) = exp{-A; } 


= ap Dinca 
= apd Z(n5s, {Gy i} 


Here, J, and Jọ are the transmitted and incident laser 
intensities and A, represents the spectral absorbance. n; is 
the number density of species i, L is the optical path 
length, a; is absorption coefficient, S; is the temperature 
dependent absorption linestrength of the absorption line /, 
and G,,;; is the line-broadening function. Voigt profile 
was used for the line-broadening function. Three absorp- 
tion lines located at 1388.135nm(#1), 1388.326nm(#2), 
and 1388.454nm (#3) were selected to measure tempera- 
ture and H,O concentration. Fig. 1(a) and Fig. 1(b) show 
theoretical H,O absorption spectra calculated using the 
HITRAN 2008 database [10] at 300K, 500K, 600K and 
1000K. Fig. 3 shows the temperature dependences of the 
three absorption lines. In order to increase the linearity 


(2) 
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Fig. 2 Relative intensity of theoretical H2O absorption spectra 
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Fig. 3 Relative intensity dependencies on temperature varia- 
tions for three absorption lines of H,O. 
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between the relative intensity and the temperature, the 
ratios between the intensity at 1388.135nm (#1) and the 
one at 1388.326nm as shown in Fig. 4 were used for the 
calculation of temperature. Thus, this allows for reducing 
the values of temperature error induced by the tomo- 
graphic algorithms. 


Experiments and Temperature Reconstructions 


The experimental data was obtained by the experi- 
mental setup [7] as shown in Fig. 5. A DFB laser (NTT 
Electronics Co., NUKIESGAAA) at 1388 nm with scan- 
ning range of 0.6nm was used to detect H,O vapor ab- 
sorption spectra. The laser wavelength was scanned at 1 
kHz and absorption spectra were measured to calculate 
instantaneous 2-dimensioonal temperature distributions 
using 8-path measurement cells as shown in Fig. 6. Laser 
beam was separated by an optical fiber splitter (OPNETI 
CO., SMF-28e 1310nm SWBC 1x16) and the separated 
laser beams were irradiated into the target area by 8 col- 
limators(STHORLABS Co., 50-1310-APC). The trans- 
mitted light intensities were detected by photodiodes 
(Hamamatsu Photonics and G8370-01), and their inten- 
sity signals were stored into the host computer. The data 
acquisition rate was set to 500 kHz (500 data points on 
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Fig. 4 Intensity ratio dependencies of the two absorption lines 
on temperature. 
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Fig.5 Experimental setup [7]. 
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Fig.6 2-dimensional cells constructed by 8-path lasers 


every | scan of absorption spectra) and temperatures in 
the measurement plane were also measured by chro- 
mel-alumel thermocouples with a diameter of 100um 
(KMT-100-100-120), in order to compare both results. 

Methane gas was used as a fuel. During all experi- 
ments, mass flow velocities of CH, and air were con- 
stantly maintained at each of 30€/min and 1020/min. The 
diameter of the burner was manufactured to 40mm, and 
the laser paths were set to 8mm spacing with each other. 
The temperature and the concentration of H2O measured 
at the burner center by sensors (thermocouple and hu- 
midity sensor) were 438K and 2%, respectively. 

Fig. 7 shows measured absorption spectra on the 8- 
path lasers across the Bunsen burner. The laser paths 3, 4, 
7 and 8 in Fig. 6 were set to high temperature regions of 
the burner. It can be seen that the ratio of line intensities 
of #1 and #2 has strong temperature dependences as al- 
ready explained in Fig. 3 and Fig. 4. The laser paths 1, 2, 
5 and 6 were set to low temperature region and the ratio 
of line intensities of #1 and #2 has lower values than 
those of the laser paths 3, 4, 7 and 8. 

These absorption values were used for temperature 
reconstructions. For temperature reconstructions, the 
MART algorithm was adopted, and its performance was 
validated. 

Absorptions of the transmitted light through the me- 
thane gas flame occur on the 8-paths. The absorption 
signals become integrated absorption values of the opti- 
cal paths [11, 12]. The integrated absorbance ‘A’ across a 
specific spectral line in the laser path j can be written as 
follows. 


A, j= Yn L jaz; (3) 


where the subscript i represent the grid point (i = 1, 2, ...), 
X is transition frequency. Because the integrated absor- 
bance is dependent on both temperature and concentra- 
tion, the temperature distribution has to be calculated by 
more than two different absorbance values. In this study, 
temperature and H2O concentration were determined at 
each analysis grids using a multifunction minimization 
method to minimize the spectral fitting error at 
1338.0-1338.6nm. This minimization method was calcu- 
lated by the following equation. 
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Fig. 7 Measured absorption spectra (H20) for 8-path lasers 
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The measured H2O absorption spectra was normalized 
and compared to those of theoretical spectra to minimize 
the MSE (Mean Squared Errors). The error was calcu- 
lated and minimized from each of all laser paths. This 
procedure is explained in Fig. 8. All values of œ, ;(i) 


Sa 


at all grids points were iteratively calculated using the 
ART and MART methods until convergence. The final 
values of Ti, j) and n(i, j) after iterations were then re- 
constructed for temperature and concentration fields. For 
data reconstructions, two algorithms, ART and MART 
were used and their performances were compared. 

Eq. (5) represents the equation used for ART (alge- 


braic reconstruction technique) algorithm. The unknown 
values of absorption coefficients (@,) are calculated im- 
plicitly. 


(k+ (k 
a (i) +H) _ a, (i) ) 
I rik 
A, j = ed a, j(i) ) Lj (5) 
I 2 
Dies 


Eq. (6) represents the equation used for MART (mul- 
tiplicative algebraic reconstruction technique) method. 
BL; 
(K+) _ (i) Avj 
CACO) 7 CA) + T ~(k) (6) 
i=l ay; (i) ; Lj 


+B 


Here, æ, ,(i) is expressed as Eq. (7). 


vl,j 

ay, Ci) =ni, /)-S;, TCG, Javi Guaj P (7) 

With initial or previous values of n(i, j) and TG, j), 
a, (i) values were iteratively calculated using Eq. (7). 


Implicit solution was performed until all values of 
æ, (i) at all grids points converged to constant values. 


k is the iteration number of MART algorithm, and £ is a 
relaxation parameter. @,; is absorption coefficient at the 
wavelength v,(1388.135nm), A,,; is the absorbance ob- 
tained from experiments. 

Iterative calculation was made for #1(v;, 1388.135nm) 
at first, next for #2(v2, 1388.326nm). When performing 
iterative calculations, the temperatures measured by 
thermocouples were used as the initial values. Then, the 
initial values of gas concentrations were obtained by the 
initial values of temperature using Eq. (3). Absorption 
coefficients were newly calculated using the absorption 
values from the horizontal laser 1 to 4 using Eq. (5). Af- 
ter iterative calculations for the horizontal lasers (1~4), 
the absorption coefficients of the vertical lasers (5~8) 
were also renewed. 

Fig. 9 shows the mean squared errors (MSE) obtained 
by ART algorithm for all iteration numbers. J values 
were changed in the range of (0.05 < f < 1). For the case 
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Fig. 8 Procedure for calculation of absorption coefficients 
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of ART algorithm, the MSE was the smallest values 
when £ was set to 0.1. The MSE was saturated within the 
iteration number 50. This implies that the iteration num- 
bers do not need to be large number if an optimal relaxa- 
tion value (£) is adopted. 

Fig. 10 shows the mean squared errors (MSE) ob- 
tained by MART algorithm for all iteration numbers. 2 
values were also changed in the range of (0.05 < £ < 1). 
For the case of ART algorithm, the MSE was the smallest 
values when J was set to 0.5. The MSE was saturated 
within the iteration number 10. This implies that the cal- 
culation time is 5 times shorter than that of ART algo- 
rithm if an optimal relaxation value (£) is adopted. 

Fig. 11 shows comparison of calculation performances 
between ART and MART algorithms. It can be inferred 
that the magnitudes of errors by MART algorithm are 
lower than those by ART algorithm with rapid calculation 
convergences. 

Fig. 12 and Fig. 13 show the distributions of tempera- 
ture and concentration obtained by MART algorithm. It 
can be seen that the peak region are biased to the edges. 
This is due to the fact that the center of flame burner 
(high temperature region) was installed at the edge loca- 
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Fig. 9 Variations of mean square error (MSE) with iteration 
numbers for relaxation parameters in ART algorithm. 
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numbers for relaxation parameters in MART algo- 
rithm. 
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Fig. 11 Comparisons between ART algorithm and MART 
algorithm 
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Fig. 13 Concentration distribution obtained by MART 
tion of the measurement apparatus as shown in Fig. 6. It 
was known that the peak values of the reconstructed 


temperature and H2O concentration were 429.7K and 
2.13%, respectively. 


Conclusions 


The MART (multiplicative algebraic reconstruction 
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technique) was newly adopted and its performances were 
validated for CI-TDLAS (computed tomography-tunable 
diode laser absorption spectroscopy). 

It was validated that the calculation errors by MART 
algorithm were less than those by ART algorithm. 

Further, the calculation speed of MART algorithm was 
5 times faster than that of ART algorithm. 

The minimum iteration number to be saturated at the 
minimum values of error was 10. This implies that sever- 
al dozens of iteration number would be enough for at- 
taining the minimum errors. 
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